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Abstract 

On the basis of the dynamical-quantization approach to open quantum systems, we 
can derive a non-Markovian Caldeira — Leggett quantum master equation as well 
as a non-Markovian quantum Smoluchowski equation in phase space. On the one 
hand, we solve our Caldeira — Leggett equation for the case of a quantum Brownian 
particle in a gravitational field. On the other hand, we solve our quantum 
Smoluchowski equation for a harmonic oscillator. In both physical situations we 
come up with the existence of a non-equilibrium thermal quantum force. Further, 
as a physical application of our quantum Smoluchowski equation we take up the 
phenomenon of escape rate of a non-inertial Brownian particle over a potential 
barrier. 
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1. Introduction 



From the mathematical point of view, the dynamics of an isolated material 
point is deterministically described by the Schrodinger equation 

dW(t) 

— = OT(t), 

where the Hamiltonian operator "K acting upon the time-dependent function <F(t) 
can carry features inherent in the particle such as its mass m, its position x or its 
momentum p. Moreover, the mathematical structure of K also relies on the 
phenomenological parameter h, dubbed the Planck constant h divided by 2n, 
which in turn is responsible for the signature of the quantum world. Most 
especially, the Schrodinger functions f(t) account for the remarkable 
phenomenon of superposition or interference of quantum states which is in the 
core of current researches on quantum computation [1]. 

Yet from the physical standpoint, a quantum system cannot be imagined as 
being in isolation from its surroundings. In truth, it only comes into being as far as 
its interaction with a certain environment (e.g., a measuring apparatus) is 
concerned [2,3]. Accordingly, a quantum system must be actually idealized as an 
open (non-Hamiltonian) quantum system comprising of a tagged particle 
immersed in a generic quantum environment. The jittering movement undergone 
by such a tagged particle is dubbed quantum Brownian motion which is to be 
mathematically described by quantum master equations of the general form 

dp{t) 

— = Xp(t), 

where the non-Hamiltonian superoperator L (the so-termed Liouvillian of the 
quantum open system) acting on the von Neumann density operator p(t) bears 
some environmental features, such as coupling constants (a sort of friction 
constant) and the fluctuation energy accounting for the existence of the quantum 
Brownian movement, as well as some properties inherent in the tagged particle 
such as its mass m, its position x, and the Planck constant h. In theory of quantum 
open systems the pivotal issue is therefore the following [4]: How can we build up 
or derive some physically meaningful non-Hamiltonian Liouville operators L? 

From the mathematical viewpoint, the superoperator L can be algebraically 
built up as Lindblad operators by beginning with a system-environment model 
Hamiltonian, and then making the Born and Markov approximations. This 
approach describes a sort of Markovian interaction between the non-Hamiltonian 
system (e.g., the Brownian particle) and its environment. In addition, it is assumed 
that the system and environment begin in a product state, i.e., they are initially 
independent and non-interacting [2 — 6]. Applications of the Lindblad formalism 
can be found, for instance, in quantum optics in which the environment is 
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represented by a quantized radiation field while the Brownian particle is deemed 
to be an atom or a molecule. The corresponding master equation is called quantum 
optical master equation [6,7]. 

Nevertheless, it has been argued that " one should not attribute fundamental 
significance to the Lindblad master equation" because "the Lindblad theory is not 
applicable in most problems of solid state physics at low temperatures for which 
neither the Born approximation is valid nor the Markov assumption holds" [8]. 
Furthermore, recent controversies [9] seem to point out the inadequacy of the 
Lindblad approach to fathoming the true physics of open quantum systems [5]. 

In a non-Lindblad perspective and on the ground of the path-integral 
approach to open quantum systems, Caldeira and Leggett [10] derived a 
Markovian master equation (the so-termed Caldeira — Leggett equation) 
describing quantum Brownian motion at high temperatures starting from a 
Hamiltonian modeling the environment as a bath of harmonic oscillators. Such a 
Caldeira — Leggett equation (CLE) has been extended by Caldeira, Cerdeira, and 
Ramaswamy [11] for any temperature but for very weak damping. Nevertheless, 
such Markovian CLEs may give rise to unphysical results, for they are not of the 
Lindblad form [3 — 6, 12]. In brief, like the Lindblad master equations, it has been 
claimed that Markovian CLEs cannot be considered as a bona fide description of 
quantum Brownian motion [8,12], albeit the high-temperature CLE has been 
employed for accounting for the decoherence process [2,13]. 

Alternatively, we have derived both Markovian Caldeira — Leggett equations 
[10,11] without alluding to the existence of an underlying Hamiltonian modeling 
the interaction process between particle and environment. Our non-Hamiltonian 
approach starts directly from the stochastic dynamics (Langevin and Fokker — 
Planck equations), thereby giving rise to quantum master equations by means of a 
quantization method we have called it dynamical quantization [5,14]. 

Thus, in the wake of our dynamical-quantization approach to open quantum 
systems, the main purpose of the present paper is to explore the physical 
significance of both non-Markovian and non-equilibrium effects upon the quantum 
Brownian motion by deriving, on the one hand, a non-Markovian Caldeira — 
Leggett equation and, on the other hand, a non-Markovian quantum Smoluchowski 
equation in the absence of inertial force. To this end, we lay out this article as 
follows: In section 2, we feature the Einstein — Langevin — Kolmogorov approach to 
open classical systems on the basis of which we can derive both the non- 
Markovian Klein — Kramers equation and the non-Markovian Smoluchowski 
equation. In section 3, we arrive at the non-Markovian Caldeira — Leggett equation 
and study the quantum Brownian motion of a particle in a gravitational field as 
well as its classical limit. Section 4 addresses the problem of quantizing the non- 
Markovian Smoluchowski equation and its classical limit, whereas section 5 takes 
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up the quantum tunneling of a non-inertial Brownian particle over a potential 
barrier. Lastly, summary and outlook are presented in section 6. In addition, three 
appendices are attached. 

2. The Einstein — Langevin — Kolmogorov approach 

Historically, Einstein [15] accounted for Brownian movement in terms of 
the time evolution of the probability distribution function (the so-called diffusion 
equation in configuration space). Then, Langevin [16] alternatively advanced a 
complementary description of Brownian motion in terms of stochastic differential 
equations (the so-termed Langevin equations [17]). From the mathematical point 
of view, Brownian motion has been looked at as a stochastic process within the 
underpinnings of theory of probability according to Kolmogorov [18]. Such 
contributions by Einstein, Langevin and Komogorov, we dub them the Einstein — 
Langevin — Kolmogorov approach to classical open systems. 

Section 2.1 considers the Langevin equation as much in the presence as in 
the absence of inertial force. Sections 2.2 and 2.3 in turn show how the non- 
Markovian Klein — Kramers and non-Markovian Smoluchowski equation can be 
derived, respectively. 

2.1. The Langevin equations 

In order to investigate the erratic motion (the so-called Brownian motion) 
of a tagged material point immersed in an arbitrary environment (a paradigmatic 
example of non-Hamiltonian open system), Langevin [16] extended the Newton's 
deterministic approach by carrying out a random generalization of the Newtonian 
dynamics through a set of stochastic differential equations (the so-called Langevin 
equations [5,17] ) 

dP dV(X) 
dX _P 

~dt~m' (2) 

The variables X = X(t) and P = P(t) are kinematical properties inherent to the 
Brownian particle whereas the variable ¥ = f(t) is intrinsic to the environment. 
The phenomenological parameter m characterizes the massive structure of the 
Brownian particle, while y and b are deemed to be positive constants responsible 
for coupling environment and Brownian particle, so that as y = b = the open 
system (1) and (2) renders isolated from its environment. 

In equation (1) the environmental force, given by F env (P,¥) = —2yP + 
bW(t), is made up by a memoryless dissipative force F d = — 2yP responsible for 
stopping the particle's motion and the so-termed Langevin force L(t) = bW{t) that 
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accounts for activating its movement. Both the phenomenological parameters 
b and y do control such an environmental influence: b is called fluctuation 
parameter while y is termed dissipation parameter. The existence of the 
environmental force F env (P, f) suggests a relationship between b and y. Further, 
we can readily check that the parameter b may be expressed in dimensions of 
[mass x length x time' 3 / 2 ], provided that <P(t) is in dimensions of [time~ 1/2 ]. 

From the mathematical viewpoint [5,17,18], the quantities X = X(i), 
P = P(t), and ¥ = f(t) in the Langevin equations (2.3) and (2.4) are viewed as 
random variables in the sense that there exists a probability distribution function 
T xpq ,(x,p,ilj, t), associated with the whole stochastic system {X,P, ¥}, expressed in 
terms of the possible values x = {x £ }, p = {p £ }, and ip = {ipi}, with i > 1, 
distributed about the sharp values q, p and cp of X, P, and V, respectively. In 
addition, in the Einstein — Langevin — Kolmogorov approach to classical open 
systems [15,16,18], the average value of any stochastic quantity A(X,P, V, t) is 
assumed to be expressed as 

(A(X,P,W,t)) = a(x,p,xfj,t)T xpy (x,p,x(j,t)dxdpdx(j, fulfilling the 
normalization condition (1) = fff*™ T X pw(x,p,xlj,t)dxdpdi(j = 1. 

As the coupling parameters do vanish, i.e., as b = y = 0, as well as assuming 
that T XP ip(x,p,xlj, t) = <5(x — q)S(p — p)8(\p — cp), the Langevin equations (1) and 
(2) become the deterministic Newton equations in terms of the average values 
(X) = q and <P> = p. 

In the absence of inertial force, i.e., as far as the condition |F env (P, f)| » 

dP 

— is concerned in (1), the focus on the motion of the Brownian particle is gained 
through the single differential equation 

dX 1 dV(X) b 

— = — H W(t) C5\ 

dt 2my dX 2my w " K J 

That is, the momentum P is said to be eliminated from the description of Brownian 

movement in the absence of inertial force. 

2.2. The non-Markovian Klein — Kramers equation 

The Langevin equations (1) and (2) give rise to the Kolmogorov equation in 
phase space [5,18 — 20] 

dT{x,p,t) 

\; - =WF{x,p,t), (4) 
ot 

where the Kolmogorovian operator K acts upon the (marginal) probability 
distribution function 7 (x,p, t) = ff™ T XPW (x,p,xp,t)dip according to 
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oo k 



KF(x,p,t) = XZ riCfe-^a^-W l^^P^Hx.p.t)] (5) 



r\(k — r)! dx k r dp r 

k=lr=0 



with 



(6) 



A^- r ' r \x,p,t) =lim 

£^>0 

In (6) both increments AX = X(t + e) - X(t) and AP = P(t + e) - P(t) are to be 
calculated from (1) and (2) from the integral equations AX = — / P(t)dt and 

Ap = -^SB. £ -2y f t+£ P(t)dt + b f t+£( P(t)dt. Moreover, the average values, 

((AX) fe ~ r )((AP) r ), are to be calculated about the sharp values q and p, i.e., 
T X p V {x,p,xp,t) = 8{x - q)S{p -p)F v (jp,t). 

The phase-space Kolmogorov equation (4) describes the time evolution of a 
Brownian particle immersed in a non-Gaussian environment since the coefficients 
(6) rely on the infinite moments of the Langevin force: (L^) ...L(t k )), with 
k = 1,2, ... , oo. According to the Pawula theorem [20], there exists no non-Gaussian 
approximation to (4) complying with the positivity of T(x,p,i). A sufficient 
condition leading to a Gaussian Kolmogorov equation (4) is then to consider 
l x 2 — x il 3 « 0- where x 2 = x(t + e) and x x = x(t) are associated with the random 
variables X(t + e) and X(i), respectively. 

So, taking into account lim f {V(t'))dt' = 0, the Brownian motion of a 

particle immersed in a generic stationary environment turns out to be described 
by the Langevin equations (1) and (2) and the corresponding Gaussian 
Kolmogorov equation (the so-called Fokker — Planck equation) 



dT _ pdT d 
dt m dx dp 



dV eff (x, t) 



+ 2yp 



d 2 T 

T + 2 Y m£(t)—, (7) 



dx 

where T = T (x, p, t) and V eii (x, t) the effective potential 

V eff (x, t) = V(x) - x^A Y mS{oo){W{t)) (8) 
with the mean value of <P(t) given by (V(t)) = lim - \l +£ {V(t'))dt' . Further, the 

£->0 £ t 

time-dependent diffusion coefficient D(t) = 2ym£ (t) is expressed in terms of the 
function 

£ W=^/(£), («) 
where /(t) is the dimensionless function 
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t + E 



/(t) = limj || (W^mt^dt^. (10) 



It is readily to check that the time-dependent function £ (t) defined by (9) has 
dimensions of energy, i.e., [mass x length 2 x time -2 ]. Hence we call it the 
diffusion energy responsible for the Brownian motion of the particle immersed in a 
generic Gaussian environment. 

The Fokker — Planck equation (7) is to be solved from a given initial 
condition T(x, p, t = 0) evolving towards a steady solution lim T(x, p, t) ~ T (x, p). 

t->oo 

In this steady regime the diffusion energy (9) displays the asymptotic behavior 

b 2 

lim^oo £ (t) « £(°°) = , provided that lim /(t) « /(°o) = 1. The physical 

4ym t— >oo 

significance of this long-time behavior has to do with the fact that environmental 
fluctuations do possess Markovian correlations, i.e., the correlational function /(t) 
exhibits a local (short) feature decaying to one in the stationary regime. By 
contrast, non-Markovian effects show up in the non-equilibrium time window 

b 2 

< t < oo. Moreover, the steady diffusion energy £ (oo) = -j— yields the 

Markovian fluctuation-dissipation relation in the form b = ^\ym£ (oo), meaning 
that the fluctuation coefficient b can be determined in terms of the friction 
constant y, the mass of the particle m, as well as £ (oo). 

The characteristic feature underlying the concept of time-dependent 
diffusion energy (9), £ (t) = £ (oo)/(t), is that it sets up a general relationship 
between fluctuation and dissipation processes as well as fulfilling the validity 
condition < £ (t) < oo. Both cases £ (t) = and £ (t) = oo should be disregarded, 
for they may violate the fluctuation — dissipation relation: the former case may 
lead to dissipation without fluctuation, while the latter one may give rise to 
fluctuation without dissipation. 

If CF(t)) = 0, /(t) = 1, and the environment is deemed to be a heat bath in 
thermodynamic equilibrium at temperature T, so that we can identify the steady 
diffusion energy £ (oo) with the thermal energy k B T, where k B is the Boltzmann 
constant (see Appendix A), then the non-Markovian Fokker — Planck equation (7) 
reduces to the Markovian Klein — Kramers equation [5,17,21]. Hence, we term (7) 
the non-Markovian Klein — Kramers equation. 



7 



2.3. The non-Markovian Smoluchowski equation 



The Langevin equation (3) gives rise to the following Kolmogorov equation 
in configuration space [5,18 — 20] 

dfix.t) v(-l) k d 



dt 

k=l 



= Z\r^ [i4fe(x ' t)/(x ' t)L 



with f(x,t)=f™ m J : xw (x,ip,t)dxp. The coefficients A k (x, t) being given by 
A k (x,t) = lim[^_>l, where the average values, ((AX) k ), are to be calculated 

e->0 L £ J 

about the sharp values q, i .e., ^F w (x, ip, t) = <5(x — q)T v (jp,t) from the Langevin 
equation (3) in the integral form: 

AX = X{t + a) - X{t) = + f§- C £ WW- Accordingly, the 

configuration-space Kolmogorov equation (11) in the Gaussian approximation, 
|x 2 — x-J 3 « 0, reduces to the following non-Markovian Fokker — Planck equation 



df{x, t) Id 



dt 2my dx 



dV ef f (x, t) 



dx 



f(x,t) 



£(oo) d z f{x,t) 



where D(t) = i s the time-dependent diffusion coefficient and V eff (x, t) the 

effective potential (8). 

For the case in which ("^(t)) = 0, /(t) = 1 as well as £ (oo) = k B T the 
Fokker — Planck equation (12) becomes the Markovian Smoluchowski equation 
[5,17,22]. For this reason, we dub (12) the non-Markovian Smoluchowski 
equation. 

We wish to point out that both the Markovian Klein — Kramers equation 
[21] and the Markovian Smoluchowski equation [22] can be obtained as a special 
case from (7) and (12), respectively, without postulating ab initio both the 
statistical properties (L(t)) = (L(t)L(t')) = 4ymk B TS(t — t') of the Langevin 
force L(t) as is commonly achieved in the literature [16,17,22]. By contrast, in our 
approach the important fact is the asymptotic behavior of (¥(t)) as well as /(t) in 
the stationary limit. In other words, both the statistical properties above turn up as 
sufficient but not necessary conditions for attaining the thermodynamic 
equilibrium state. 

Having showed how non-Markovian effects can turn up in both Klein — 
Kramers and Smoluchowski equations, in the next sections 3 and 4 we intend to 
examine how quantum effects upon these kinds of non-Markovian Fokker — Planck 
equations can arise via a non-Hamiltonian quantization process. 
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3. Deriving the non-Markovian Caldeira — Leggett equation 



Having derived the Kolmogorov equation in phase space (4) for a Brownian 
particle immersed in a non-Gaussian environment, we now wish to quantize it by 
means of a non-Hamiltonian quantization process called dynamical quantization 
[5,14]. First, we obtain the equation of motion 



dxQc, r\, t) 
dt 



= | Ky(x,p,t)e £ P"dp, (13) 



after performing upon (4) the Fourier transform j(x, n, t) = S™ m T(x, p, t)e ipv dp, 

where the exponential e ipv is deemed to be a dimensionless term. Then, the 
Kolmogorov stochastic dynamics (4) is said to be quantized by introducing into 
equation (13) the following quantization conditions 

x 1 =x + — ; x 2 = x - —, (14) 

whereby the transformation parameter h from the change of variables 
(x, r\) h> (x 1 ,x 2 ) is dubbed Planck's constant having dimensions of angular 
momentum, i.e., [mass x length 2 x time' 1 ]. The geometric meaning of the 
quantization conditions (14) lies at the existence of a minimal distance between 
the points x x and x 2 on account of the quantum nature of space, i.e., |x 2 — x x \ = 
\t]h.\, such that in the classical limit h -> 0, physically interpreted as |r;h| « 
|x 2 — x-J, the result x 2 = x x = x can be readily recovered. 

Thus, after inserting both quantization conditions (14) into (13) we arrive 
at the non-Gaussian quantum master equation 

oo 

dp(x 1 ,X 2 ,t) C fX 1 -X 2 x fp (*i-*2) 



dt 



= J Ky(^-y^,p,t)e tp ^^dp, (15) 



which describes the quantum Brownian motion of a particle in the presence of a 
generic quantum fluid, for instance. 

On going from the classical equation of motion (13) to the quantum 
dynamics (15) via quantization conditions (14), we have replaced the classical 
function x = xi x > V> with the quantum function p = p(x 1; x 2 , t), since it turns out 
to depend on the Planck constant, h. We dub p(x 1( x 2 , t) the von Neumann function. 

Taking the Gaussianity condition |x 2 — xj 3 «0 into account, our non- 
Gaussian quantum master equation (15) reduces to 
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(x 1 -x 2 ) 2 p, (16) 



dx 7 / 



where the effective potentials V e ^?(x fc , t), with k = 1,2, are given by V^(x k , t) = 

^( x /c) — x k j4ym£ h (oo)(W(t)) and DjiCO i s the time-dependent quantum diffusion 
coefficient T> h (t) = 2ym£ h (oo)/(t). 

It is relevant to notice that in the Gaussian quantum master equation (16), 
the time evolution parameter t, the mass m, the frictional constant y, as well as the 
both functions CF(t)) and /(t) are deemed to be non-quantized quantities, that is, 
they are h-independent, whereas the classical diffusion energy £ (oo) has been 
subject to a quantization process, i.e., £ (oo) -> £ h (oo). 

The master equation (16) describes the quantum Brownian motion of a 
particle immersed in a generic Gaussian quantum environment. If, for the specific 
case of a heat bath consisting of quantum harmonic oscillators with oscillation 
frequency o>, we could identify the Brownian particle's steady diffusion energy 
£ h (oo) with the mean thermal energy of the bath in thermodynamic equilibrium, 

then we obtain £fr(°°) = — — = — coth (see Appendix A). Accordingly, the 

equilibrium quantum fluctuation — dissipation relation reads 

b h = I lymoSh coth reducing to b = ^/4ym£(oo) in the classical limit, i.e., as 

k B T » o>h/2. Moreover, the quantum diffusion constant D h (t) in (16) reads 
2?h(0 = ywz&jhcoth (——)/(0' yielding at high temperatures, T » a)h/2k B , the 

\2kgT/ 

classical diffusion constant in (7) for thermal systems: D(t) = 2ymk B TI(t). On the 
other hand, the zero-point diffusion constant reads T)^~°\co) = ymoSh. Thus, for 
thermal quantum systems, the non-Markovian quantum master equation (16) 
reads 

* Tt = K ff Vi. o - v e y (* 2 , t) j P _ _ _ — j 

(dp dp\ ( o>h \ 

-ihyfe - x 2 ) ^ - - iym*> coth {—) /(Ofe -x 2 ) p (17) 

which describes a Brownian particle subject to averaging effects present in the 
effective potentials v!S(jx k , t) as well as non-Markovian effects through the 
correlational function /(t). In addition, it is worth underscoring that (17) has been 
derived for any initial condition p(x 1; x 2 , t = 0). 

As long as CF(t)) = and /(t) = 1, our quantum master equation (17) 
leads to the Markovian Caldeira — Leggett equation 
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fa—=[V(x 1 ,t)-V{x 2 ,t)]p-, 

(dp dp\ 
ihy(x 1 — x 2 ) I — I — iymoy coth 



2m \dxf dxl ) 



h 2 fd 2 p _d 2 p\ 



—i 




found out by Caldeira, Cerdeira, and Ramaswamy [11] following the Feynman path 
integral formalism and making assumptions on the weakness of the damping 
Y « of - Yet our derivation has shown that Feynman's formalism as well as such an 
assumption are thoroughly unnecessary features for reaching (18). Moreover, on 
the condition that the thermal reservoir holds at high temperatures, i.e., 
coth(a>h/2/c B r) ~2/c B r/a>h, our master equation (17), along with (¥(t)) = and 
/(t) = 1, yields the following master equation 



originally found by Caldeira and Leggett [10] upon making use of the path integral 
techniques and assuming that particle and environment are initially uncorrelated. 
This assumption is non-realistic and leads to non-physical results [12,23]. Our 
approach on the contrary has shown that the high-temperature Markovian 
Caldeira — Leggett equation (19) can be derived for any initial condition as long as 
averaging and non-Markovian effects could be neglected in our master equation 
(17) at high temperatures. 

In addition, it has been claimed that the Markovian Caldeira — Leggett 
equation (19) may also lead to unphysical results because it is not of the Lindblad 
form [12]. To overcome this difficulty, terms has been added to (19) for fitting with 
the Lindblad framework. Yet, according to our approach such an ad hoc procedure 
is unnecessary given that our quantum master equation (17) conveys non- 
Markovian effects, so bypassing Lindblad's condition. 

Lastly, because our non-Markovian quantum master equation (17) contains 
both the Markovian Caldeira — Leggett equations (18) and (19) as special cases, we 
dub it the non-Markovian Caldeira — Leggett equation. 

3.1. A quantum Brownian particle in a gravitational field 

In order to provide physical significance to our non-Markovian Caldeira — 
Leggett equation (17), let us consider the correlational function of the form 
/(t) = 1 — e _t,/tc (see Appendix B) and the case of the quantum Brownian motion 
of a particle of mass m moving due to the gravitational energy, V(x) = mgx, where 
g denotes the value of gravity acceleration. This free -falling quantum Brownian 
motion is described by the following dynamics in quantum phase space 





(19) 
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dW 



1 

m 



mg 
2y 



WO) 



_ + 27 _ [pW]+Dh(0 __ 



(20) 



,m£h(00) <«?(*)> and the time- 



in terms of the effective momentum p = p' H 

dependent quantum diffusion coefficient V h (t) = 2ym£ ^(oo) ^1 — e ^j. The 
function W = VK(x, p, t) in (20) is the Wigner transform 



oo 

If r?h r?h 

x,p,t)=- j p(x + y ,x- y , 



(21) 



performed upon the von Neuman function in the master equation (17). 

To solve equation (20) we start with the following initial condition that 
couples Brownian particle with the environment 



1 Jl 

W(x,p,t = 0) = — e V' 
nn 



r p 2 mx 2 
mh fir, 



) 



(22) 



where x r = (4y) 1 denotes a sort of relaxation time. It is straightforward to check 
that (22) complies with the Heisenberg fluctuation relation AP(0)AX(0) = h/2. 

Taking W(p,t) = f™ oo W(x,p,t)dx into account, equation (20) turns out to 
be rewritten as the quantum Rayleigh equation 



dW(p,t) d d 2 W(p,t) 
the time-dependent solution of which reads 



(23) 



W{p,t) = 



e 4 ^(0 



with 



(24) 



mh _1 mcoh / o)h \r __1 t c ( -— -— W 
\W=-e . + _coth(— ^ + ^^(. *)].<2S) 

The probability distribution function (24) leads to <P> = and (P 2 ) = 2A h (t). 
Hence, the momentum fluctuation AP(t) = j(P 2 ) — (P) 2 reads 



mh __1 mcoh / wh \r __1 t c / -— -— \1 
AP(t) = J- e , r+ _coth(— - +J —^(e <r-e *)].(26) 
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Because AP(t) = AP'(t) we conclude that the mean value of V as well as the 
gravity effect present in the quantum Klein — Kramers equation (20) are 
unobservable. Further, it is readily to check that non-Markovian effects account for 
diminishing the strength of the fluctuation (26). 

From (24) in the steady regime we can derive 

p 2 



W{p) = 1 — e mMh (A), (27) 

jnma)h. coth \^~fj 

meaning that both relaxation and non-Markovian effects in (24) die out in thermal 
equilibrium. In other words, both relaxation and non-Markovian features in (24) 
are intrinsically non-equilibrium effects. 

Furthermore, the root mean square momentum (26) gives rise to the 
following thermal quantum force 

dAP(t) , 

X(t) = = -2yV2yhmB(t) (28) 

with the dimensionless function B(t) 
B(t) 

/ coh \ / — — \ — 

t 2 co coth j ) y e Tr — e t cj + (t c — r r )e T r 



(t c — T r ) |(t c — T r )e T r + T r O) COth ( 2^7 ) [ T r ^ e %r ~ 1 j + — 



(29) 

displaying the following physically accessible time scales: the evolution time t (the 
observation time, for instance), the correlation time t c , the relaxation time r r , the 
oscillation time t osc = co -1 , as well as the quantum time t q = h/k B T. If B(t) > 0, 
then the force (28) is said to be attractive: K(i) < 0. Otherwise, if B(t) < 0, then 
(28) renders repulsive, i.e., K(i) > 0. 

The quantum force (28) is a non-equilibrium effect since it vanishes in the 
steady regime, whereas at short times t « t c , t r the force (28) becomes attractive 
with the constant value K(t) = —ly^lyhm. This fact implies that the root mean 
square momentum (26) is a differentiable quantity for all time t > 0. 

By way of example, we take |B(t)|~l, reckon with the Munro and Gardiner's 
values for the parameters y, h, and m in Ref. [12], i.e., y~10 11 s" 1 , 
h~10" 34 m 2 kgs _1 , and m~l(T 26 kg, and then obtain the magnitude |5C(t)|~10~ 13 N. 
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This simple numerical example suggests that the strength of the quantum thermal 
force (28) could be measured in experiences, for instance, using trapped ions [24] 
in which measurement of forces of order of yoctonewton, i.e., 1CT 24 N, has been 
recently reported. 

At zero temperature, i.e., coth(a)h/2k B T) ~1, the dimensionless factor (29) 
becomes 



B^=°)(t) = 



(e tr-e ^) + (t c -T r ) 



t_ 

e T r 



:,(30) 



(t c — r r ) |(t c — r r )e T r + a)T r |r r T r — lj + t c (l — e t c 



while at high temperatures, i.e., as coth(coh/2k B T) ~2k B T/a)h, it turns out to be 



2k, 



B(t) = 



,Tt? (e T r - e ^ + h(t c - r r ) 



t_ 

e T r 



(31) 



h(t c - T r ) jh(t c - r r )e T r + 2k B Tx r [r r (e T r - + t c - e f c 



thereby suggesting that the non-equilibrium quantum effect (28) could show up at 
both high and zero temperatures. 



3.2. Classical limit 

The classical limit h -» should be physically interpreted as the classical 
thermal energy is deemed to be too large in comparison to the quantal energy, i.e., 

k B T » so that £ji(°o) = ~f~ cot h ~£(°°) = k B T. So, from the quantum 

Wigner function (24) in the steady regime at high temperatures, we can derive the 
well-known Maxwell — Boltzmann probability distribution function /"(p) = 



, 1 e 2mk B T for the effective momentum p = p' + — — \ mkBT (y(t)) while in 

j2nmk B T 2y v y 



the classical limit the sort of quantum force (28) becomes 

F(t) = -2Y^Jmk B TG(t) 
with the dimensionless function G(t) being 

t r f r — e t^j 



(32) 



G(t) = 



(33) 



(t c -t r ) 



^ei-ljl^l-ei)] 
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Taking into account |G(t)|~l, y-lO^s -1 , m~10~ 26 kg, /c B ~10~ 23 m 2 kgs" 2 K -1 , and 
r~1000K, we find |F(t)|~10 _22 N, which is the magnitude of the classical thermal 
force (32) exerted by a heat bath at 1000K upon the Brownian motion of a free- 
falling particle of mass 10~ 26 kg . 

The dimensionless function (33) approximates to G(t)~ — (2yt c ) _1/2 at 
short times. Hence the non-equilibrium classical thermal force (32) reads 

F(t ) = prnkBT blQws 

up in the Markovian limit t c -> 0. This upshot 

reveals the pivotal importance of non-Markovian effects for the existence of the 
physical concept of non-equilibrium thermal force (32) in the classical realm. In 
mathematical parlance, non-Markovian features account for the differentiability 

property of the root mean square momentum AP(t) = yj{P 2 ) — (P) 2 for all time t 
in the classical domain. Accordingly, the Langevin equation 

dp i , ^ 

— = -2yP + j4ymk B TV(t), 

underlying the non-Markovian classical Rayleigh equation 

df(p,t) d d 2 f(p,t) 

obtained as classical limit of (23), should be interpreted as a differential stochastic 
equation and not as an integral equation according to a determined rule of 
interpretation (Doob's interpretation, for instance. See Coffey et al.'s book in Ref. 
[17] as well as [25]). 

3.3. Discussions 

The predominant paradigm for describing quantum Brownian motion can be 
summarized as follows [5,6,8]: Envisage an environment consisting of a set of 
harmonic oscillators coupled to a Brownian particle; having built up the 
Hamiltonian or Lagrangian function of the whole system (particle plus 
environment), canonical quantization procedure (Dirac or Heisenberg 
quantization) could be employed to enter into quantum world. Once established 
this quantum Hamiltonian picture portraying dissipation along with fluctuation, 
Feynman's path integral formalism or projection operator techniques may be then 
invocated for deriving master equations for the density matrix describing the 
quantum motion of the Brownian particle, after getting rid of the medium's 
variables. On the basis of the former procedure, Caldeira and Leggett [10] 
quantized the Markovian Klein — Kramers equation [21] and obtained the so-called 
Markovian Caldeira — Leggett equation. Yet, this quantum master equation is 
plagued with the problem of positivity for it is not of the Lindblad form [2]. 
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By contrast, according to our dynamical-quantization approach to quantum 
Brownian motion we have derived the non-Gaussian Kolmogorov quantum master 
equation (15) without building up a Hamiltonian model for the isolated system 
(Brownian particle plus environment). Our quantization method is intrinsically 
non-Hamiltonian in full. By the same token, we reckon that our Caldeira — Leggett 
quantum master equation (17) eschews the issue of complying with the Lindblad 
requirement because it gets the non-Markovian gist of the interaction between the 
Brownian particle and the thermal reservoir in the Gaussian approximation. 

4. The non-Markovian quantum Smoluchowski equation 

In order to provide some generality to our method of non-Hamiltonian 
quantization, we turn to quantize the non-Markovian Smoluchowski equation (12) 
describing a harmonic oscillator. The most general case given by the quantization 
of the Kolmogorov equation in configuration space (11) is taken up in Appendix C. 

4.1. The harmonic oscillator 

We start with (12) for V(x) = kx 2 /2 at points x x and x 2 , i.e., 



d*(*i. _. _J d _ 

dt 2my dx 1 

dx(x 2 ,t) 1 d 



dVeffC*!, t) 



dt 



2my dx 2 



dx x 
dV en (x 2 ,t) 



dx 7 



X(xi> 

x(.x 2 , 



_ d 2 y(x 1; t) 
-, \ d 2 x(x 2 ,t) 



kxt 

where V eK (x k , t), with k = 1,2, are the effective potentials V eK (x k , t) 
2x k ^2myS{'X)){ l ¥{t)), whereas D(t) is the time-dependent diffusion coefficient 
D(t) = 2D(t) = "^7^(0 associated with both solutions x( x i> = Jf( x i> an d 
X(x 2 , t) = yjf{x 2 ,t). 

Multiplying (34) and (35) by j(x 2 ,t) and x(. x i>t)> respectively, and then 
adding the resulting equations we arrive at 



dt 



—\— 

2my \px x 



dV eK (x lt t) 



dx 1 



d 


dV ef{ (x 2 ,t) 


dx 2 


dx 2 



+ 2)(t) 



d 2 % d 2 % 
+ 



dx 1 



dxn 2 



(36) 



where f = ^ r (x 1 ,x 2 , t) = x(. x i> x(. x 2> 0- By quantizing via (14) in the form 
x x = x' + y and x 2 = x' — y, as well as making use of the following relations 



a _ i a i a a 

a*! 2 h dr\ dx 2 



-— + -—, the classical equation of motion (36) 

2 dxi hdrj n v J 



changes into the quantum equation of motion 



dp k 

-£ = — P + 

dt my 



kx dp 



+ 





£ h (0 


d 2 p 


dri 


2my 


dx 2 



+ 



h 2 dt] 2 



(37) 
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expressed in terms of the effective position x = x' — ^2my£ li (oo)(<f / (t)). 

Upon going from (36) to (37) we have replaced the solution ^ = %(x lt x 2 ,t) 
with the h-dependent function p = p(x, r\, t) and employed the subscript notation 

to display the quantum nature of the diffusion energy 

which turns out to be now expressed in terms of the Planck constant h. Notice that 

__t_ 

both the constants k and y as well as the dimensionless function lit) = 1 — e *c 
hold classical during the classical-to-quantum transition. 

By making use of the Fourier transform 

oo 

W(x, V ,t) = — j pix.rj.^e^drj, (38) 

— oo 

which changes the variables from quantum configuration space (x,vjh) onto 
quantum phase space (x, p; h), equation of motion (37) turns out to be written 
down as 



dW k dW k dW £ h (t)d 2 W 

= X V 1 h 

dt 2my dx 2my dp 2my dx 2 



2£ h (t) 2 
p z 



2my myh 2 



W (39) 



in terms of the phase-space function W = W(x, p, t). Because the exponential 
factor e ipv in (38) is to be a dimensionless term, it follows that the variable p is to 
have dimensions of linear momentum. 

4. 1. 1. The initial condition 

Our quantum phase-space Smoluchowski equation (39) may be solved 
starting from the non-thermal initial condition 

1 ( a P 2 i x2 \ 

W{x,p,t = 0) = — e V h ^) (40) 
nn 

leading to the distribution F(x, p) = <5(x)<5(p) in the classical limit, h -> 0. The 
constant a has dimensions of time per mass, that is, we may set a = 1/my, for 
instance. The Gaussian probability distribution function (40) gives rise to both 

variances (P 2 (0)) = ^ and (X 2 (0))=— , which fulfill the Heisenberg 
fluctuation principle V(P 2 (0))(X 2 (0)) = \. 

4.1.2. The time-dependent solution 

A time-dependent solution to (39) reads 

W{x,p,t) = — e~t^ Bh(t) e"ra (41) 
nn 
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with 



BM =^ e ^-H 1 -' tr+ T^rX e tr ~ e '<)} (42) 

expressed in terms of the evolution time t, the relaxation time t r = my/k, as well 
as the correlation time t c . 

Solution (41) yields <P(t)> = (X(t)) = and the following fluctuations 

AP(t) = V(P 2 (t))-(P(t)) 2 = -JL=, (43) 



AX(t)=V^ 2 (t))-(X(t))2= jM> f (44) 



satisfying the Heisenberg constraint AP(t)AX(t) = \ The fluctuation (43) gives 



rise to the concept of thermal force 



dAP(t) h d5 h (t) 

nt) = -^-= (45) 

[2S h (t)]2 



whereas (44) generates the thermal velocity 



dAX(t) 1 dB h (t) 
Kt) = — ^ = r — £A (46) 

at 2[2S h (t)]z ar 



where 



d5 h (t) fch -L 2£ h (a>) t, 



dt m 2 7 2 my t 



t c ( -± 



e tr+ — — e tc-e ^ J. (47) 



h _ , dB h (t) 



my dt 



feh 



Since S h (t = 0) = — and — f^- = — , both quantities (45) and (46) are 



t=o m 2 y 2 ' 



well defined at t = 0, \.e.,T(t = 0) = - —and V(t = 0) = -— —.This 

v y 2 -\J 2my v J 2my yj 2my 

upshot implies that both fluctuations (43) and (44) are differentiable functions at 
t = 0. 

4. 1.3. The long-time regime 

In the long-time regime t -> oo, physically interpreted as t » t c , t r , the 
probability distribution function (41) renders steady, i.e, 

2£ h (°°)p 2 



Ws t (x,p)=— e h 2 e 2f h(°°), (48) 
nn 
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h I k 

while the fluctuations (43) and (44) become, respectively, AP(oo) = - I— — and 

2 \| £h(°°) 



AZ(oo) = J^^. In contrast, we find out J'(oo) = V(oo) = 0, meaning that both 
quantities (45) and (46) are non-equilibrium effects. 

4.2. Thermal systems 

4.2.1. Quantum limit 

We now assume that our Brownian harmonic oscillator is immersed in a 
heat bath consisting of quantum harmonic oscillators having the same oscillation 
frequency co, so that the Brownian particle's quantum diffusion energy £fr(°°) can 

be identified with the heat bath's quantum thermal energy E = — coth (-7^-), i.e., 

£ h (oo) = — coth (-7^-) [see Appendix A]. In the quantum limit, i.e., at zero 
temperature, the thermal diffusion energy £h(°°) should be replaced with 
£ h (oo) = ^ in the formulae above. 

4.2.2. Classical limit 

In the classical limit h -» 0, the time-dependent solution (41) becomes 
T(x, p, t) oc S(p)e fi w that in turn leads to the function 

f(x,t) = f T{x,p,t)dp = e~W), (49) 

—00 v 

which is a solution of the non-Markovian Smoluchowski equation (12) for classical 
thermal systems. In (49), S(t) is the function 



B(0 = 2k ° T 



_t_ 1 / _t_ __L\"i 
1-e tr+ c _ (e f r- e t c \ I (50) 



obtained from (42) as lim h ^ S h (t) = S(t). 

In the classical limit the root mean square momentum (43) vanishes, 
whereas (46) becomes 



AX(t) = 



k B T 



t 



tr 



1 c £ > 4 \ (e - e (51) 



breaking down the Heisenberg constraint, i.e., AP(t)AX(t) = 0, as expected. This 
implies that in the classical realm the momentum becomes a deterministic variable 
while the displacement holds as a stochastic one. In other words, the momentum 
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variable is said to be eliminated from the classical description of Brownian motion 
in the absence of inertial force. 

It is worth stressing that the root mean square displacement (51) is 

— - — . Only in the Markovian limit t c -> 0, 



2myt c 



(51) renders nondifferentiable. 



Furthermore, we can derive the Maxwell — Boltzmann distribution 

kx 2 

1 2 r 

function /(x) = -j==e 2k s T for the thermal position x = x' — - A /2m7/c B r( ( F(t)) 
as the classical limit of the (marginal) probability distribution function (48). 



5. Application: quantum tunneling 

On account of the environmental fluctuations, the escape rate of a Brownian 
particle over a barrier separating two metastable states — in a double-well 
potential, for instance — is known as the Kramers problem [17,21,26]. In the 
classical domain, the transition of such an particle over the potential barrier from a 
metastable state toward another state is said to be an activation phenomenon, 
while the corresponding metastability phenomenon in the quantum domain 
characterizes a tunneling process. In both realms of physics, escape rate 
phenomena are non-equilibrium effects taking place in the (Markovian) steady 
regime. 

In 1940 Kramers [21] extended the Einstein — Langevin — Kolmogorov 
approach to study the issue of metastability of a Brownian particle in phase space 
and found an escape rate characterized by non-quantum, thermal, and steady 
effects. In [27] we have generalized Kramers' theory to open classical systems far 
from thermal equilibrium. In this section, we wish to extend Kramers' technique to 
our quantum Smoluhowski equation (39) reckoning with two physical situations: 
Section 5.1 deals with non-thermal systems and section 5.2 considers thermal 
systems. Section 5.3 brings some discussions. 

5.1. Non-thermal open quantum systems 

In order to calculate the quantum Kramers escape rate of a non-inertial 
Brownian particle over a potential barrier, we consider as a starting point the 
steady solution (48) in the form 

1 \ IM , 2£ h (oo)p 2 l 

W st (x,p) =—e ** 2 J. (52) 
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We perform an expansion of the potential function V(x) = kx 2 /2 in a Taylor series 
around a certain point x , i.e., V(x)~V(x Q ) + (/c/2)(x — x ) 2 , and then construct 
the steady non-equilibrium function 



W ne Jx,p) = a<p(x,p)e 



[H*o)+§Q-*o) 2 ] 2£ h (oo) 



h 2 k 



(53) 



a being a constant and <p(x, p) a function to be determined. On inserting the 
function (53) into our quantum Smoluchowski equation (39), given by 



k(x - x ) dW neq kp dW neq _ f h (oo) d 2 W r 



'neq 

dx 



2 dp 



+ 



dx< 



+ 



k 2£ h (oo) 



h 2 



w neq = o, 



(54) 



with W neq = VK neq (x, p), we obtain the equation of motion 



"neq 



d 2 (p{x,p) d(p{x,p) d(p{x,p) 
£h(°°) — — 2 k(x - Xq) — — kp — — — = 0, 



dx 

which changes into 



dx 



dp 



(55) 



d<p($) 



in terms of the new variable £ given by 

f = /?(x-x ) - p. 



(56) 



(57) 



The constant /? in (56) and (57) is to have dimensions of mass per time. Hence, we 
may put /? = my/2, for instance. A solution to the differential equation (56) reads 



-k 



; zke 
J 



(58) 



27r£ h (oo) 

— oo 

fulfilling the condition <p(f -> oo) = 1. This result (58) is only possible if the 
potential curvature is negative, i.e., k < 0. Hence, we hereafter use x = x b , and 
k = —ma) b , where the quantity a) b denotes the particle's oscillation frequency over 
the potential barrier at x b . So the barrier top is located at point x b , while the two 
bottom wells are at x a and x c , such that V(x c ) = = V(x a ), with x a <x b . 



Using (58), the stationary function (53) becomes 



2aa) h 



£ h (oo) 



Y A /27rmf h (oo) 
from which we can find out the probability current as 



2£h(oo) 2 
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oo 

=-\ 

m J 



P^neqfo = X b ,p) dp = - AC — ,__ == = , (60) 



V(x b ) 

where we used the result 

oo f=p 

1 / 7T vl/2 



— oo 



it and v being non-negative constants. 

The number of particles through the well located at x a is calculated as 
n a = JJ.qo W4q( x <P)dpdx = arch, whereby we used the equilibrium distribution 

4£ h (oo)p 2 mo)| ^ 2 

function VK eq (x,p) = ere h2m&) a 4£ h(°°) a coming from the non-equilibrum 
function (53) for <p(_x, p) = 1 calculated at the well x a , with k a = ma)a. 

Making use of (60) and n a = anh, the quantum Kramers escape rate of a 
Brownian particle immersed in a non-thermal environment is then written down 
as 

r(oo) = — = - — — — - (61) 

n a 47r£ h (oo)^ h2<u 4 - Y 2 £^(po) 

fulfilling the condition < y£fr(°o) < hco 2 ,. 

It is worth highlighting that our main upshot (61) does display dimensions 
of inverse of time, as it should be expected, as well as being compatible with the 
quantum dissipation-fluctuation relation £fr(°°)- 

Our quantum escape rate (61) decays exponentially with respect to the 
ratio— V(x b )/£ h (co), In addition, the so-called quantum Arrhenius factor is given 
b = _J: <h 2 

5.2. Thermal open quantum systems 



For thermal systems, i.e., £fr(°°) = — coth (-^—), our quantum Kramers rate (61) 

2 \2knT J 



<2k B TJ 

becomes 



2V(x b ) 



w fr hC0th (5?r) 
- coth (^)>^ 2 -r 2coth2 (^) 



O) 2 

r O) = 7T7^-T= == (62) 
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which is valid for low temperatures: ycoth(-^— ) < 2a) b . At zero temperature 

\2/CgT/ 

equation (62) reads 

2V(x b ) 

r v=*K<*) = ^-^=, (63) 

with 7 < 2(1);,. Both results (62) and (63) show how dissipative effects affect the 
quantum tunneling process in the Smoluchowski regime. This dissipative tunneling 
phenomenon is in accordance with Ankerhold's complaints [28]. It is worth 
noticing that dissipative effects in both (62) and (63) account for enhancing the 
quantum tunneling rate. This outcome is in contrast to the classical Kramers 
escape rate [21] 

a> a a> b 

T(oo) = e k B T 

Any 

in which damping effects are responsible for diminishing the Arrhenius factor. 
5.3. Discussions 

Because of the wide practical and theoretical interest in chemical physics, 
for instance, recent studies have been come out showing how quantum effects 
upon the Smoluchowski equation and the Kramers rate can arise. 

First, on the basis of the Caldeira — Leggett Hamiltonian system-plus- 
reservoir model [5,8,10] in tandem with the path integral formalism, Ankerhold 
and co-workers [8,29] have derived a quantum Smoluchowski equation (QSE) and 
explored its physical meaning applying it to chemical reactions, mesoscopic 
physics, and charge transfer in molecules, for instance. Meanwhile, further works 
[30] have revealed some drawbacks underlying such a QSE, for it may violate the 
Second Law of Thermodynamics, for example. It is relevant to point out that 
Ankerhold's QSE for the case of a free Brownian particle coincides with the 
classical Smoluchowski equation . 

Another Hamiltonian system-plus-reservoir account aiming at to quantize 
the Smoluchowski equation has been developed in [31]. Here a QSE, which for free 
Brownian motion differs from the classical Smoluchowski equation, is derived and 
the problem of quantum tunneling is studied at zero temperature, in contrast to 
Ankerhold's survey. This Indian group's approach does not depend on the path 
integral formalism but is based on the canonical quantization. 

Coffey and co-workers [32] in turn have obtained a QSE in the semi-classical 
region of quantum mechanics on the basis of Brinkman's hierarchy in terms of the 
time evolution for the Wigner function obtained from the Caldeira — Leggett 
Hamiltonian system-plus-reservoir model. Although Coffey's QSE can be applied to 
the dynamics of a point Josephson junction, for instance, it is plagued with the 
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same disease as Ankerhold's QSE, namely, both quantum and classical dynamics 
for free Brownian motion are identical. 

Starting from the quantum linear Boltzmann equation, Vacchini [33] has 
derived a QSE in which the dynamics of a free quantum particle is distinct from a 
classical free particle. 

Lastly, Tsekov [34] has arrived at a version of QSE without referring to any 
quantization process. He starts with a Hamiltonian system-plus-reservoir model 
described by a non-linear Schrodinger equation in the Madelung representation of 
quantum mechanics and then obtains a QSE for a free Brownian particle different 
from the classical Smoluchowski equation, too. 

Alternatively, in previous works [5,14] and in the present study we have 
derived a quantum Smoluchowski equation based on a non-Hamiltonian 
quantization method without making use of the canonical quantization as well as 
the path integral formalism. Our QSE (39) for a free particle (k = 0) does not 
coincides with the classical Smoluchowski equation. Moreover, our equation (39) 
can be applied to looking at quantum tunneling process at zero temperature (see 
result (63)). 



6. Summary and outlook 

In this paper we have taken up the dynamical-quantization approach to 
quantum open systems without assuming the existence of a Hamiltonian 
framework underlying the interaction between a Brownian particle and its 
environment. Our non-Hamiltonian approach features the following upshots: 

(a) By quantizing the non-Gaussian Kolmogorov equation in phase space 
(4), we have derived the quantum master equation (15) describing the quantum 
Brownian motion of a particle immersed in a generic environment (e.g., a non- 
thermal quantum fluid). For the case of a thermal quantum reservoir, our master 
equation (15) in the Gaussian approximation yields the non-Markovian Caldeira — 
Leggett equation (17) that has been solved for a particle in a gravitational field. 
Here, we have put forward the concept of thermal quantum force (28) that is a 
non-equilibrium effect carrying both non-Markovian and relaxation features. 
Moreover, we have shown that Markovian effects account for the differentiability 
property of the root mean square momentum (26) in the classical limit. 

(b) Quantizing the non-Gaussian Kolmogorov equation in configuration 
space (11) leads to the quantum equation of motion (C.4) reducing to (C.6)in the 
Gaussian approximation (see Appendix C). For the case of a harmonic oscillator 
immersed in a general non-thermal environment, our Gaussian quantum equation 
(C.6) gives rise to the non-Markovian quantum Smoluchowski equation in phase 
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space (39) that has been solved. For a thermal quantum heat bath we have derived 
the thermal quantum force (45) as well as the thermal quantum velocity (46) that 
are well defined at t > and at zero temperature. In the classical limit (in the high- 
temperature regime) the root mean square momentum (43) vanishes meaning 
that the momentum variable is a deterministic one in the non-inertial classical 
Brownian motion. Moreover, the classical limit of (43) yields the quantity (51) that 
is differentiable for all time owing to non-Markovian effects. By contrast, in the 
Markovian limit the root mean square displacement (51) renders non- 
differentiable, as it is well known since Einstein's 1905 paper. 

(c) Our quantum Smoluchowski equation (39) has been also applied to 
describing the tunneling process of a non-inertial quantum Brownian particle in a 
non-thermal dissipative medium. Our main finding is the quantum Kramers rate 
(61). For thermal reservoir (61) becomes (62) that in turn is valid in the low- 
temperature regime (including the zero temperature case). According to our 
quantum escape rates (62) and (63) dissipative effects account for enhancing the 
quantum tunneling rate. 

Lastly, in a forthcoming article we intend to tackle the problem of quantum 
Brownian motion in the presence of bosonic and fermionic heat baths as well as 
non-Gaussian environments as described by our quantum Kolmogorov equations 
(15) and (C.4). 
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Appendix A. The physis of the environment 

Imagine that the environment is made up of a set of quantum harmonic 
oscillators having the same oscillation frequency o>. It is readily to show that the 
mean energy of this quantum heat bath after attaining the thermodynamic 
equilibrium at temperature T is given by [35] 



where k B is dubbed Boltzmann's constant. The h-dependent energy, coh/2, 
corresponds to the zero point energy of the heat bath at zero temperature and 
k B T denotes the classical thermal energy of the heat bath at high temperature, i.e., 
T » a)h/2k B . 
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Appendix B. The non-Markovian correlational function 

The correlational function /(t) is given by 

t+E 

/(t) = lim- [[ (V(t')V(t"))dt'dt" . 

£^0 £ J J 
t 

Under the condition lim /(t) « 7(cx>) = 1, the autocorrelation function 

t->oo 

(W(t')W(t")) can be built up as 

/ (t/+t/Q \ 

(nt'W)) = U-e 2tc )ff(t'-t"), 

where t c is the correlation time of ¥(£) at times t' and t". Accordingly, we find the 
non-Markovian correlational function 

_t_ 

/(t) = 1 - e tc, 

reducing to /(t) = 1 in the Markovian limit t c -> 0. 
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Appendix C. Quantizing the Kolmogorov equation in configuration space 

We wish to quantize the non-Gaussian Kolmogorov equation in configuration 
space (11). To this end, let Xi = x(. x i> be a solution of the non-Gaussian 
Kolmogorov equation (11) and j 2 = j(x 2 , t) another one atx 2 : 



-1^ 



dt k\ dx k 

k=l 1 



(-l) fc d k r - , 

ji\A k (x 1 ,t)x(x 1 ,t)\, 



(CI) 



dt 



=T-=i 



(-i) fc d k r _ , 

k [A k (x 2 ,t)x(x 2 ,t)\. 



i=i 



dx" 



(C.2) 



where A k (x lt t) and A k (x 2 , t) are coefficients associated with solutions j(x 1( t) and 
j(x 2 ,t), respectively. Multiplying (C.l) and (C.2) by j(x 2 ,t) and j(x 1( t), 
respectively, and then adding the resulting equations we arrive at 



d$(x 1 ,x 2 ,i) 
dt 



k=l 



(-l) fe ( d k - d k - \ 

nAy-ix^t) + -^ k A k (x 2 ,t) )f(x lt x 2 ,t) 



dx k 



(C.3) 



where $(x 1 ,x 2 ,t) = j(x 1( t) j(x 2 , t) = V/(*2, t)/(*i,0- B Y quantizing via (14) 

d id Id d id id 

and making use of the relations - — = -- — -—and - — = -— + -—, 

b dx-t 2dx hdri dx 2 2 dx h dt} 

the classical equation (C.3) becomes the quantum Kolmogorov equation in 
configuration space 



dp(x,r],t) 
dt 



= L h p(x,r],t) 



(C.4) 



with the h-dependent operator L h given by 



U - V \( _A_ _ ilX A 00 ( x _^l t) + (— + --) A™ (x + V - t] 

Lh ~L k\ \2dx hd V ) Ak V 2 ,t ) + {2dx + hdr ] ) Ak \ X+ 2 ' l ) 

k=l 

(C.5) 

The quantum Kolmorov equation (C.4) exhibits non-Gaussian features in full. Yet, 
in the Gaussian approximation |x 2 — xj 3 « 0, it changes into 

dp 1 d 2 V^(x,t) 1 dV^{x,t)dp r, d 2 V^{x,t)dp 

-p + - =^ h 



dt my 



dx 2 

£h(°°) 
my 



/(t) 



2my dx dx 2my 



d 2 p 4 d 2 p 



dx< 



d V 



+ 



dx 2 h 2 dry 



(C.6) 



with the effective potential 
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V™(x, t) = V(x) - x^AymS h {oo){W{t)) (C. 7) 

Considering (C.6) at point x' for a harmonic potential V(x') = kx' 2 /2, it reduces to 
equation (37) in terms of the effective position x = x' — ^2my£ h( 00 )( ( ^(t)). 
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